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We consider the conductance of a quantum tube connected to a metallic contact. 
The number of angular momentum states that the tube can support depends on 
the strength of the radial confinement. Wc calculate the transmission coefficients 
which yield the conductance via the Landaucr formula, and discuss the relation of 
our results to armchair carbon nanotubes embedded in a metal. For Al and Au 
contacts and tubes with a realistic radial confinement we find that the transmission 
can be close to unity corresponding to a contact resistance close to h/2e 2 per band 
at the Fermi level in the carbon nanotube. 

© 1999 Academic Press Limited 
1. Introduction 

Since the recent discovery of the carbon nanotubes by Iijima there has been a significant 
progress in the studies of the conducting properties of both single-walled || and multi-walled 
[Q carbon nanotubes. Conductance of a mesoscopic system connected to metallic reservoirs is well 
understood and is usually described by the Landauer formula [§]. For quantum point contacts in 
semiconductor structures and in metallic nanowires it is well establish experimentally that the 
differential conductance to a good approximation is quantized in units of Go = 2e 2 /h and at zero 
temperature given by G = GqN where N is the number of propagating modes. In carbon nanotubes 
with metallic contacts most experiments show that the conductance is less than the conductance 
which one should expect for a smooth interface between tube and metal, e.g. G — 4e 2 /h for metallic 
single-walled tubes where the extra factor of 2 comes from the two ir bands that are crossing the 
Fermi level 0. The reasons for this lower conductance are still not fully known. Theoretically, 
several groups have considered the effects of vacancies [Q, disorder distortion and doping 
pp[ on the conductance of carbon nanotubes. The conducting properties have also been studied 
in e.g. the context of junctions between different metallic carbon nanotubes jO], Aharonow-Bohm 
effect in the presence of a magnetic field |l2j , and the Luttinger liquid behavior of a one-dimensional 
gas of interacting electrons p3|. Also the ideal "hollow quantum cylinder", i.e. a two-dimensional 
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electron gas on a cylinder, has been studied in context of the difference between strip-like wires and 
tubes |Q. However, with the exception of the recent qualitative study of Tersoff [^5j and the recent 
modeling- works of Anantram et al. Jl6| and Sanvito et al. |l7| |, less attention has been focused on 
the conditions for a good transmission between tube and a metal contact which is an important 
issue for practical devices with carbon nanotubes, or other quantum tubes. 

In quantum point contacts an adiabatic interface between the wire and reservoirs ensures a trans- 
mission coefficient close to unity [fl8) . The condition for adiabaticity is that the shape of the contact 
region varies slowly on the scale of the Fermi wave length. In the opposite case with an abrupt 
interface, i.e. quasi-one-dimensional lead connected to a wide two-dimensional contact, Szafer and 
Stone |]l9f found that the transmission rapidly increases to unity as the width of the confined region 
exceeds half of the Fermi wavelength, thus giving a reflectionless contact. 

For the contact between a quantum tube and a three-dimensional metal it is not obvious that the 
assumption of an ideal reflectionless contact applies and the aim of this work is to study the contact 
resistance for this case. 

The model we are studying is that of a hollow quantum cylinder of radius R T contacted by a three- 
dimensional free-electron metal which we for convenience model by a cylindrical wire with radius 
R c 3> R T , see Fig. [l| The system thus has full cylindrical symmetry and the angular momentum 
quantum number m can therefore be used to label the scattering states. 

For the coupling of the quantum tube to the contact it is necessary to take a radial confinement 
potential for the quantum tube into account and here we model the confinement by an attractive 
delta-function potential. As an example we apply this model to metal contacts of Al or Au; the 
quantum tube parameters are chosen to mimic armchair carbon nanotubes: the strength of the 
confinement can be related to the work function for the material that constitutes the tube and in 
the case of a carbon nanotube we relate it to the work function of graphene. It should be noted 
that the employed free electron model does not fully describe the actual band structure of carbon 
nanotubes. Nevertheless, a study of contact resistance within this idealized model should yield 
valuable insights which are relevant to real materials. 

The paper is organized as follows: In Section II the eigenstates of a quantum tube connected 
to a cylindrical metal contact are found. In Section III these eigenstates are used to construct the 
scattering states to find the transmission coefficient, and hence the conductance of the contact. In 
Section IV, we apply our model to contacts between an armchair carbon nanotube and a metal. 
Finally, in Section V discussion and conclusions are given. Essential details of analytical calculations 
are given in Appendices A and B. 

2. The eigenstates 

We separate the discussion into two parts: first we find the eigenstates in the tubular geometry 
and then the eigenstates for the cylindrical metal contact. In Section III the matching of these 
eigenstates are used to construct the scattering states of the contact. 

2.1. Quantum tube 



The quantum tube of radius R T with otherwise free electrons is modeled by the Hamiltonian 
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with a confining potential given by an attractive delta function potential 

Vr(r) = -H8(r - Rr), (2) 

where the confinement strength H is taken positive. 

The eigenstates of the Schrodinger equation have the form 

% m (r,(f),z) = R m (r)Xm(<P)ip m (z), (3) 

with angular and longitudinal wave functions 

Xm (<P) = (27r)- 1 /2 exp[zm0]j (4) 

Tpm{z) = [k m (E)}- 1/2 exp[±ik m (E)z], (5) 

1/2 

where the angular momentum quantum numbers m are integers, k m (E) = [ 2 jp L {E — e m )] is the 
wave vector associated to the longitudinal free propagation, and E = E m +e m is the total energy of 
the state. Here, E m > is the energy associated to the longitudinal propagation and e m < is the 
(binding) energy associated to the transverse motion. We can relate the strength of the confinement 
to the work function W = \e m \ — fi 2 k 2 ?/2m e: which is the energy required to remove an electron at 

— 1/2 

the Fermi level (disregarding surface charge effects). The normalization [k m (Ej] is chosen such 
that the propagating modes carry the same amount of current. 
The radial wave function R m (r) satisfies 



2 d d 

T h r — 

dr 2 dr 



2rrt e e r , 
U 2 



Rm(r) = -jRr6(r - R T )R m (r), (6) 



where 7 = 2m e HR T /h 2 is a dimensionless confinement strength. For the bound states (e m < 0) and 
r 7^ R T this equation has the form of Bessel's modified differential equation p0[ . The solutions are 
given by modified Bessel functions of order in of the first and second kind, so that the full solution 
is given by 

td ( \ J ^m-^m(^m^) : f ^ Rt /^\ 

m[r) ~ 1 B m K m {K m r) , t > Rt [ ' 

1/2 

where n m = \2m e \e m \ /Ti\ . At r — R T , the radial wave function is continuous and the appropriate 
matching condition for the derivative dR m (r) /dr at r — R T is found by integrating Eq. (|^) from 
i?~ = R T — + to i?+ = R T + + . In this way the matching conditions become 



R m (R^) — R m (R T ) — 0, 



dR m (r) 



dr 



dR m (r) 



dr 



= -7- 



R m (R T ) 
R T 



and we get the following equation for the normalization coefficients 

-^m(^m-^T) ^■m\^"mRr) 

K m R- 



An— l(^m-^r) ~t~ {^tuRt) k ^ r> ^ Km— l(^m-^T) ~t~ -^m-fl (^m-^T 



(8) 
(9) 






(10) 



Non-trivial solutions exist if the determinant vanishes, and hereby the wave vector n m is a solution 
to the equation 



(11) 
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where the result for the Wronskian W {K m (x) ;I m (x)} = l/x has been used (^TJ. Expanding Eq. 
( p"T| ) in the small-K m i? T limit pc[ we find that a bound state with angular momentum mh exists for 
7 > 2m. The number of bound states for a certain value of 7 is given by N — Int(7/2) + 1 where 
Int(x) is the integer part of x. Thus, there is always at least a single bound state corresponding to 
m = 0. 

From the matching conditions and the normalization of R m (see Appendix [A|) , it follows that 

ImjKmRr) TS / \ p , , (12) 

K m (K m R T ) ^m{K m r) , K T < r 



V2 



Am 1 R 



I m (f^rn Rt: 



-1/2 



(13) 



A plot of the radial wave function is shown in the inset of Fig. 3. Increasing the confinement strength, 
the radial wave function becomes more localized which is accompanied by an increase in the binding 
energy. 

2.2. Metal contacts 

For the metal it is convenient to assume a cylindrical geometry and consider a cylindrical wire of 
radius R c 3> i? T - The Hamiltonian is written as 



2m e 

with a hard-wall confining potential 



d 2 d 2 Id Id 2 
dz 2 dr 2 r dr r 2 dd> 2 



+ V c (r), (14) 



V c (r) = { V ° ' T D <Rc . (15) 
w \ oo , R c < r K ' 

Obviously, for r > R c , \P(r, <f>, z) = and for r < R c the eigenstates have the form 

V um (r,<p,z) = R vm (r)Xm(<P)ipum(z), (16) 

Ruin \J) — CurnJm ( Kv m r), (17) 

Xm {4>) = (2vr)- 1 / 2 exp[ l m ( /»], (18) 
^m(z) = [k l/m (E)]~ 1/2 exp[±ik vm (E)z], (19) 
where J m is a Bessel function of the first kind of order to, n 2 m ~ 2m e £ mn /h 2 is a wave vector 

1 /2 

corresponding to the radial energy e um , k um {E) = \^f- (E + V£ — e um )\ is the wave vector of 
the longitudinal motion, and E is the total energy of the state. Again the normalization [k um (E)] 1 
makes the propagating modes carry the same amount of current. 

The boundary condition for the radial wave function leads to J m {ni>mRc) = 0, from which we 
find k„ numerically. Since J m (x) ~ (2/ttx) 1 ^ 2 cos(x — mn/2 — n/4) for large x Q, we have 
HumRc — [y + m/2 — l/4)7r with v = 1, 2, 3, . . .. The normalization C vm is given by (see Appendix 

0) 



C vm = — / - , (20) 

Rc\ Jm~l{ K vmRc)Jm+l{K'iimRc) 



which is a real number. 
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3. Transmission of contact 

We consider an electron in the tube in mode to incident on the contact (see Fig. [l]) and compute 
the transmission and reflection coefficients. We construct the scattering states in the basis of the 
eigenstates of the Schrodinger equation (see previous section). Since the angular momentum is a 
conserved quantity, the transmitted and reflected parts of the wave function also have the same 
quantum number to. In the quantum tube (z < 0) the scattering state is given by 



* m (r,0,z) = 
and in the contact [z > 0) by 



Rm{r) 



exp(ifc m z) 



+ r r 



exp(— ik n 



i(r, <f>, z) = t vm R um (r) 



exp(im0) exp(ifc„ TO 2r) 

V 27T 



(21) 



(22) 



Here r m is the reflection amplitude for mode to and t vm is the corresponding transmission amplitude. 
We assume that the effective electron mass is the same in the two materials so that the continuity of 
^m(f, 4>, z) and d^mi^r, 0, z) jdz at z — are appropriate boundary conditions. For carbon nanotubes 
and metals like Al and Au this is a reasonable approximation (to* ~ 1.2 too pl| ). For general details 
on how to account for differences in the effective mass and the underlying symmetry of the lattice 
we refer to Refs. and references therein. The boundary conditions lead to 



v— 1 Qvm 
is i n 



2ft 



2—iv=l "nr. 



(23) 
(24) 



where g um = \J k um /k m (R m \R vm ) , with the radial overlap defined as (R m \Rum) — 
J °° &rrRm{r)R vm {r). In addition we have the sum-rule Y^=i {Rm \Rum) 2 = 1, which can be used 
to verify the numerical convergence. The overlap can be calculated analytically (see Appendix [^) 
and the squared overlap is given by 



(Rm \Rur. 



[^m(^i'm^r) ^vmlm^mR^Km^rnRc^Jm+l^vmRc^ 



(K m i? T ) + (Kj/ m i? T ) ['An-l( K i'm-Rc)>4i+l( K i/in^c)] 



(25) 



[Kmi^mR^Im—l (/%n-^r)-^m+l (^)71^t) I^i {^mRT^jKm — 1 {j^mR r r)^m+\ (^m^r)] 

The total transmission from mode to in the quantum tube into the contact is thus 



T 

+ in 



^ y l^vm \tv 

2 



1 + E^U oir. 



(26) 



where V V m projects onto the propagating modes (k vm real) of the metal contact. Here we have 
assumed that the lengths of the quantum tube and the contact are semi-infinite so that tunneling 



6 



Superlattices and Micro structures, Vol. ??, No. ??, 1999 



through evanescent modes can be neglected. These should be included in the case of two metal 
contacts connected by a quantum tube of finite length. Introducing real and imaginary parts by 
E^li Qlm = T ',n + we obtain 

(i + r^) 2 + r^ 2 

The reflection probability lZ m can be calculated in a similar manner which provides us with the 
usual sum-rule T m + lZ m = 1, ensuring the conservation of probability current density. 

To summarize, the transmission probability of mode m can be calculated from Eq. ( p7f ) with 

k m (E) = [^(E- e m )] 1/2 and k vm {E) = (E + V° - e vm )] 1/2 , with e m = ^K 2 m being the 
energy of the mth transverse mode in the tube and similarly e vm — ^-K^ m is the transverse energy 
in the contact. Here n m and k„ are solutions to I m (K m R T )K m (n m R T ) = 7 _1 and J m (K>umRc) = 0, 
respectively. For a numerical implementation, an upper cut-off v c in the sum over modes in the 
contact is needed and the sum-rule for the squared radial overlap is then a measure of the numerical 
convergence for a given cut-off. 

When choosing values for the confinement parameters, 7 and V®, we take into account that the 
Fermi momenta of the quantum tube and the metal can be different since the relevant electrons 
are located at the Fermi levels of the two materials. For the metal contact we use the known Fermi 
energies for e.g. Al and Au to relate the confinement potential to the Fermi level as Ep = E + V£, 
with the Fermi energy being defined positive. For the tube, the Fermi energy enters as Ep = E+\e m \. 
When the two materials are brought into contact the chemical potentials align, but the difference in 
Fermi wave vectors remains. Thus for the metal contact E + V® = h 2 [fcp] 2 /2m e , and for the tube 
E+\e m \ = h 2 [fc£] 2 /2m e . For the tube we need to specify 7, which follows from the work function 
W = \e m \ — h 2 [fcp] 2 /2m e . We have neglected the charge density induced at the interface by a 
mismatch of the work functions. For a discussion of this in the context of the screening properties 
of one-dimensional systems, see e.g. Ref. E3]. 



4. Contact resistance of single-walled armchair carbon nanotubes 

The (n, n) armchair single- walled carbon nanotube can be regarded as the result of rolling one 
sheet of graphite (with the carbon atoms in a hexagonal lattice) in the direction of one of the bonds 
[pl| . The resulting tube has a periodicity a ~ 0.246 nm along the tube axis (z-axis) and a radius 
R T ~ n x ^/3a/2iT with An atoms along the perimeter, arranged in two rows that resemble a chain 
of armchairs, see Fig. ||. Their metallic character is caused by two ir bands crossing the Fermi level 
at a wave vector hp ~ 2n/3a. 

As discussed recently by Tersoff Jll| the metallic armchair carbon nanotubes have electrons at 
the Fermi level which can be regarded as having an angular momentum quantum number m = 0. 
In order to apply our simple model to the problem of the contact resistance of (n, n) single-walled 
carbon nanotubes embedded in a free-electron metal we notice that kpR T ~ n/\/3. In Fig. || the 
transmission probability at the Fermi level is shown for several values of kpR T corresponding to 
(n, n) armchair carbon nanotubes for various values of the dimensionless confinement strength 7. In 
the particular case of an Al contact, the mismatch is given by fcp/fcp 1 ~ 0.49 and the corresponding 
transmission is presented in panel (a) of Fig. In panel (b) we show similar results for an Au 
contact for which fcp"/fcp u ~ 0.70. 

In order to estimate 7 we relate it to the work function of the carbon nanotube which is of the 
order 4 — 5 eV. For the quantum tube we associate a work function to the m = bound state via 
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its binding energy, i.e. W = -^r- — Ep where k,q is the solution to Io(kqR t )Kq(kqR t ) = r y^ 1 . In 
Fig. ^ this work function is shown as a function of the confinement strength for quantum tubes 
corresponding to (n, n) armchair carbon nanotubes. From this we estimate that 7 < 10 — 20 is 
a reasonable regime for the curves shown in Figs. [| and [|. This means that a transmission close 
to unity (per band at the Fermi level) can be expected if the armchair nanotube is embedded in 
free-electron metals like Al or Au. Assuming a work function W = 4.5 eV of the nanotube we have 
used the corresponding value of 7 to calculate the transmission for the different nanotubes. For 
(n, n) armchair nanotubes with Al and Au contacts we find T Kl ~ 0.93 and T An ~ 0.98, respectively. 
In the case of matching Fermi wave vectors (fcp/fcp = 1) we get T ~ 0.87. For the considered tubes 
(3 < n < 17), these transmissions are found to be almost independent of the specific value of the 
tube indices (n, n). 

Here we have only taken the geometry- related contact scattering into account. Physically, a lower 
transmission can be caused by electrons being scattered by interface imperfections/roughness, de- 
viations from a spherical Fermi surface of the metal contact, and scattering due to non-matching 
work functions of the nanotube and metal. Also scattering due to the non-matching Fermi veloc- 
ities of the nanotube and the metal could be expected. However, as shown in Fig. ^, a mismatch 
between Fermi wave vectors can actually in some cases increase the transmission (and thereby the 
conductance) due to quantum interferences, even though the mismatch by itself is known to give 
rise to momentum relaxation and thereby resistance. 



5. Discussion and conclusion 

We have considered the contact resistance (in terms of transmission) of a quantum tube embedded 
in a free-electron metal. For the quantum tube we have modeled the radial confinement of the 
electron motion by an attractive delta function potential which gives rise to at least one bound 
state in the radial direction. The strength of the attractive potential can phenomenologically be 
associated to the work function of the quantum tube. Within this model we have calculated the 
transmission of a quantum tube contacted by a free-electron metal. Due to the cylindrical geometry 
of the contact, considerable analytical progress was possible and with the resulting equations the 
scattering problem is readily solved numerically. 

As an application we have considered the transparency of contacts with armchair carbon nan- 
otubes embedded in free-electron metals. Our calculations show that in the absence of scattering 
mechanisms associated to e.g. interface imperfections/roughness, deviations from a spherical Fermi 
surface of the metal contact, and scattering due to non-matching work functions of the nanotube and 
metal, the geometry itself allows for a high transparent contact between armchair carbon nanotubes 
and free-electron metal contacts. Furthermore, from this simple model we find that Al would be a 
good candidate for such a metal as it was suggested recently by Tersoff |ll|. For Au however, we 
find that the present 3D geometry allows for good contact in contrast to Tersoff 's findings for Au, 
which were based on ID considerations. 
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A. Normalization of radial wave functions 
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From the radial wave function of the quantum tube, Eq. ([T2|), it follows that the normalization is 
given by 



Am — 



K m H T 7-2 / n \ 



da aK m (a) 



-1/2 



and since |24| 



daal^(a) = a 2 [l m (a) - J m _i(a)7 m+ i(a)] /2, 
daa<(a) = a 2 [Jf* (a) - K ro _ 1 (a)iT m+1 (a)] /2, 



(28) 

(29) 
(30) 



we get the result in Eq. (jl3|). Similarly, from the radial wave function of the free-electron metal 
contact, Eq. (0), it follows that the normalization is given by 



Rc 



1 -1/2 



drr.J m (n um r) 



and since Jm(K vm R c ) = and 



daaJ m (a) = a 2 [J m {a) - J m _x(a) J m+1 (a)] /2, 



(31) 



(32) 



we obtain Eq. ( P0|) 



B. Overlap of radial wave functions 

The overlap of radial wave functions can be written as 



A n 



dv" Tl m {K m T)J m {K um T S ) H — - , / df"!'Km( K m r )Jm( K iym' r ) 

-Km(K m rt T J Ji{ T 



^mCym 7 m (K^ m i? T ) -(- KumRcIm {^mRr^-^m (^m-^cj^m+l (^vm-Rc) 



where we have used the integrals ]24 

drr/ m (ar)J m (/3r) = 
dr 7\ff m (ar) J m ((3r) = 



r {a/ m+ 1 (or) J m (/3r) + /?/ m (or) J m+ i (/?r ) } 
a 2 + /3 2 

r {aAT m+ i (ar)J m (/3r) - (3K m (ar) J m+1 {(3r)} 



(33) 



a 2 +/? 2 



(34) 
(35) 



together with the boundary condition i?„ m (i? c ) = 0. 
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Fig. 1. Contact between a quantum tube (z < 0) of radius R T and a three-dimensional cylindrical free-electron 
metallic wire (z > 0) of radius Rc S> -Rt- 




Fig. 2. A sheet of graphite (upper panel) which can be rolled up to a (n, n) armchair single- walled carbon nanotubc 
(lower panel). The tube has a periodicity a ~ 0.246 nm along the tube axis (z-axis) and a radius R T — n X \/3a/2n 
with 4n carbon atoms along the perimeter, arranged in two rows that resemble a chain of armchairs. The shown 
example is a (5, 5) armchair nanotube. 
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Fig. 3. Transmission probability T from an (n,n) armchair carbon nanotube of radius i? T — n/\/3fcp into a 
cylindrical contact with radius Rc 2> Rt- T is shown for several values of the confinement strength 7. The results 
marked with Cg> are the specific values of R T kp = n/\/3 corresponding to (n,n) tubes and the dashed lines are 
calculated curves which are shown as guides to the eye. The inset shows the radial wave function of the quantum 
tube, Eq. (fLSj), for m = 0. The nanotube and the metal are assumed to have the same Fermi wave vectors, so that 

^f/^f ~ ^' 
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Fig. 4. The same calculation as in Fig.bl but, with a mismatch between the Fermi wave vectors of the carbon 
nanotube and the metal. Panel (a) is for fcp/fcp 1 ~ 0.49 corresponding to nanotubes embedded in an Al contact and 
panel (b) is for fcp/fcp u ~ 0.70 corresponding to nanotubes embedded in an Au contact. 
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Fig. 5. Work function W of quantum tubes corresponding to (n, n) armchair carbon nanotube as a function of 
the dimensionless confinement strength 7. Experimentally, the work function of a carbon nanotube is W ~ 4 — 5 eV 
which for the shown nanotubes means that 7 < 10 — 20 is the important parameter range. 



